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FOREWORD 


The  final  technical  report  and  this  Supplement  were 
prepared  by  F.M.  White,  R.C.  Lessmann,  and  G.H.  Christoph 
of  the  Department  of  Mechanical  Engineering  and  Applied 
Mechanics  of  the  University  of  Rhode  Island  under  Contract 
F33615-71-C-158E ,  "Analysis  of  the  Turbulent  Boundary  Layer 
in  Axisymmetric  and  Three-Dimensional  Flows." 

The  contract  was  initiated  under  Project  No.  1426,  "Ex¬ 
perimental  Simulation  of  Flight  Mechanics,"  Task  No.  142604, 
"Theory  of  Dynamic  Simulation  of  Flight  Environment."  The 
work  was  administered  by  the  Air  Force  Flight  Dynamics  Labor¬ 
atory.  Wright-Patterson  Air  Force  Base,  Ohio,  Dr.  Tames  T. 

Van  Keren  fFX) ,  Project  Engineer. 

'The  work  was  accomplished  during  the  period  1  June  1971 
through  30  June  1972. 

The  report  was  submitted  by  the  authors  in  July  1972. 

This  Supplement  contains  an  outline  of  the  theoretical 
method  for  engineering  use. 

This  Supplement  has  been  reviewed  and  is  approved. 
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Diaecsions . 


It  is  thv  purpose  of  this  report  to  outline  the  coapiete  details 
for  engineering  use  of  White's  tvo-diaensional  theory  or  turbulent 
skin  friction  under  arbitrary  coapressible  flow  conditions.  This 
new  aetbod  is  not  only  the  simplest  but  also  the  aost  accurate  com¬ 
putational  scheme  in  the  literature.  This  procedure  is  recoaaended 
for  general  use  by  engineering  designers.  The  aethod  concerns  skin 
friction  only.  This  report  presents  methods  for  hand  or  digital 
I  computer  computation.  Several  ezaaples  are  presented. 
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I  INTRODUCTION' 


The  authors  have  been  engaged  for  several  years  in  the  development 
of  a  new  method  of  analysis  of  turbulent  skin  friction  under  fairlv 
arbitrary  conditions  of  freestream  and  wall  parameters.  Early 
efforts  [1,2,3]  were  devoted  to  two-dimensional  skin  friction 
calculations  .while  later  parsers  consider  axisvmmetric  skin  friction 
tto,5],  three-dimensional  skin  friction  [5],  wall  heat  transfer  [5], 
and  wall  roughness  and  transpiration  conditions  [6].  It  is  the 
purpose  of  this  report  to  outline  the  complete  details  for  engineering 
use  of  our  two-dimensional  theory  [3]  of  turbulent  skin  friction 
under  arbitrary  compressible  flow  conditions.  It  appears  to  the 
authors  that  this  new  method  is  not  only  the  simplest  hut  also  the 
most  accurate  computational  scheme  in  the  literature.  Therefore  we 
feel  justified  in  recommending  this  procedure  for  peneral  use  bv 
engineering  designers  - 

This  methoc  concerns  skin  friction  only.  Ko  attempt  is  made  to 
estimate  the  value  of  anv  "traditional"  integral  parameter  such  as 
a  shape  factor  or  a  momentum  thickness.  Such  parameters  could  easilv 
be  added  to  the  present  analysis,  but  we  adamantly  discourage  their 
use.  In  our  opinion,  shape  factors  and  integral  thicknesses  are 
of  little  significance  and  serve  mainly  to  divert  the  attention  of 
the  engineer  from  his  basic  problem,  the  turbulent  skin  friction. 

This  th^me,  admittedly  biased,  dominates  all  our  work  [1-6]. 
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li  '"‘•I  T I XI;  OF  THE  BASIC  THEORY 


Ir.  reference  [3],  the  use  of  a  compressible  law-of-the-wall 
c^nd  a  Crocco  temperature  approximation  are  shown  to  lead  to  closure 
of  the  boundary  layer  continuity  and  momentum  eauations  for  two- 
dimensional  compressible  turbulent  flow.  The  variables  which  arise 
quite  naturally  from  this  analysis  are  the  skin  friction  variable: 

X  =  <2/C,)1/2  ,  C  =  2r  /o  U2  ,  (1) 

f  f  wee 

which  varies  with  the  dimensionless  coordinate  in  the  freestreais 
direction: 

x*  =  x/L  (2) 

A  basic  parameter  in  the  theory  is  the  dimensionless  freestream 
velocity  distribution: 

V  =  U  (x)/U  =  V(x*)  (3) 

e  o 

The  quantities  L  and  Uq  are  reference  constants;  they  may  be  taken 
equal  to  unity  if  one  is  prepared  to  handle  the  re  'ulting  "unit  Reynolds 
numbers”  and  "unit  lengths"  in  the  basic  equatici.  Ir  ^;tr  <■>•-•*1 
we  commonly  take  L  equal  to  the  body  length  ho  that  x*  varies  from 
zero  to  unity)  and  take  Uq  equal  to  the  initial  velocity  lMx=0). 

We  will  illustrate  both  these  approaches  in  what  follows,  sine.,  it 
is  our  experience  that  improper  use  of  "x*"  and  in  th?  theory 

is  the  main  cause  of  erroneous  results  when  the  theory  is  used. 

? 


The  analysis  not  only  requires  accurate  knowledge  of  V(x*)  but 
also  its  first  and  second  derivatives,  V*  and  V".  Therefore  it  is 
essential  that  an  accurate  curve-fit  or  other  smooth  formula  be 
provided  for  the  distribution  V(x*).  Failure  to  provide  a  smooth 
curve-fit  is  the  second  most  common  source  of  error  in  the  theory. 

We  shall  also  illustrate  this  problem  with  examples.  Note  that  there 
are  no  limitations  upon  tMx),  which  may  be  subsonic,  transonic, 
supersonic,  or  hypersonic.  The  theory  is  valid  for  all. 

The  second  important  runnirg  parameter  in  the  theory  is  a  sort 
of  "stretched”  Reynolds  number,  I*: 

R*  =  Rl/(1/V)'  ,  (4) 

l  /? 

where  JL  =  (U  L/v  )(u  /v  )(T  ft  )  . 

X  o  e  e  w  e  w 

The  prime  indicates  differentiation  with  respect  to  x*.  The  Revnolds 

numoer  is  a  constant  for  low  speed  adiabatic  flow  [ref.  1]  but 

may  vary  with  x*  at  high  speeds  due  to  wall  and  freestream  temperature 

variations.  Finally,  the  freestream  Mach  number  variation  Mg(x*)  and 

the  wall  temperature  ratio  T  (x)/T  (x)  are  combined  into  a  sir.eie 

w  e 

parameter  A(:<*),  first  used  by  van  Driest  [7]  _n  a  flat  plate  analysis: 

A  =  /Ur  7t  -  l)/[sin  ^(a/c)  ♦  sin  ^"(b/c)]  , 

3W  C 


(5) 


where  T  is  the  adiabatic  wall  temperature  and  a/c,  b/c  vary  oi.^v 
aw  " 

with  Mach  number  and  temperature  ratio: 

a  =  ?  +  T  -  2T  ;  b  =  T  -  T  ;  c2  =  (T  +T  )2  -  4T  T  .  (K) 

aw  we  aww  aw  v  we 

Numerical  values  of  this  parameter  A  are  shown  in  Figure  1-  This 

figure  may  be  used  to  pick  off  values  of  A  when  making  a  hand  computation 

with  the  theory.  Otherwise,  when  using  the  comouter  program  in  this 

report,  computation  of  A  from  '5)  and  (6)  is  already  programmed. 

Note  that  A  varies  only  slowly  with  Mach  number  ar.d  temperature 

ratio.  Therefore,  an  accurate  curve-fit  expression  for  M  (x*>  and 

e 

(Tw/Te)  is  not  needed;  any  reasonable  approximation  will  do. 

All  of  these  parameters  arise  in  the  basic  differential  equation 
for  computing  A(x*)  from  an  arbitrary  (known)  distribution  of  V,  M  , 

and  T  /T  .  The  ecuation  is: 

w  e 


dA 

dx* 


a/VV  ti  ,  9  a-2  r*0.07]  +  (1 /vr;[3  a2  r,  p*0*0”7] 


Tim 


0.16  f*  A' 


(7) 


For  low  speed  flow,  A  =  1.0  and  £a.(7)  reduces  tc  the  incompressible 
analysis  of  reference  [1].  The  functions  f*  and  g*  arise  from 
the  integral  coefficients  of  [1,3]  and  varv  only  with  the  ouantitv 
(A/A  ).  Thev  are  plotted  in  Figure  2,  which  is  the  second  figure 

TTV3.X 

required  for  making  a  hand  computation  from  the  theory.  Finally,  the 

auantitv  "A  "  is  the  value  of  A  at  which  boundary  iaver  separation 
max 
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Figure  1.  The  Van  Driest  Parameter  A,  from  Equations  (5,6) 


gure  2.  The  Functions  f*  and  g*  for  Use  with  Equation  (7). 
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occurs.  This  is  correlated  in  [1,3]  with  R*  and  A  ,  as  follows: 


\nax  =  8'7A  <0' 

As  this  point, f*  approaches  zero  (Figure  2)  and  thus  the  derivative 

2 

(dX/dx*)  becomes  infinite  and  =  -/X  approaches  zero.  This  of 

course  is  the  definition  of  the  separation  point.  The  theory  at 

present  cannot  be  used  beyond  the  separation  point. 

Numerical  values  of  f*  and  g* ,  taken  from  [1],  are  given  in 

Table  1.  In  the  computer  program  which  follows,  those  values  are 

curve-fit  with  the  following  expressions  which  give  excellent  accuracy 

in  the  range  0.36  <  X/X  <  1.0: 

**  max  ■" 

f*  =  (2.434  Z  +  1.443  7?)  exp[-  44.0  Z6]  ; 

(9) 

g*  =  1  -  2.3  Z  +  1.76  Z3 

where  Z  =  1  -  (X/X  ). 

max 

If  (X/X  )  is  less  than  0.36,  then  f*  =  g*  =  0,  and  Eq.(7)  is  not 
max 

necessary.  The  same  is  true  if  R*  is  negative,  which  corresponds  to 

a  favorable  pressure  gradient.  In  these  cases,  the  following  simplified 

differential  equation  is  used: 

(  X  <  0.36  X  or  R*  <0) : 
max 


TABLE  1 


THE  FUNCTIONS  f*  AND  ft* 
FOP  USE  WITH  EQUATION  (7) 


A/A 

max 

g* 

f* 

<  0.36: 

Use  Equation  (10)  instead. 

0.36 

0.00237 

0.07575 

0.38 

0.00436 

0.13522 

0.40 

0.00774 

0.22410 

0.42 

0.01310 

0.34015 

0.44 

0.02098 

0.47212 

0.46 

0.03176 

0.60479 

0.48 

0.04562 

0.72521 

0.50 

0.06255 

0.82546 

0.52 

0.08242 

0.90229 

0.54 

0.10503 

0.95566 

0.56 

0.13018 

0.98735 

0.58 

0.15763 

0.99997 

0.60 

0.18715 

0.99642 

0.62 

0.21853 

0.97944 

0.64 

0.25159 

0.95155 

0.66 

0.28615 

0.91498 

0.68 

0.32204 

0.87164 

0.70 

0.35914 

0.82313 

0.72 

0.39732 

0.77079 

0.74 

0.43646 

0.71574 

0.76 

0.47646 

0.65887 

0.78 

0.51725 

0.60092 

0.80 

0.55873 

0.54250 

0.82 

0.60085 

0.48407 

0.84 

0.64353 

0.42601 

0.86 

0.68673 

0.36864 

0.88 

0.73040 

0.31219 

0.90 

0.77449 

0.25682 

0.92 

0.81895 

0.20268 

0.94 

0.86377 

0.14986 

0.96 

0.90890 

0.09844 

0.98 

0.95431 

0.04846 

1.00 

1.00000 

0.00000 
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Ill  HAND  CALCULATION* 


The  basic  relation,  Eq . ( 7) ,  is  quite  suitable  for  hand  calculation 
and  its  integration  may  be  acconclished  either  numerically  or 
graphically.  The  following  steps  outline  a  tynical  portion  of 
the  calculation: 

1)  given  tne  value  of  X  at  some  intermediate  position  x*. 

2)  Establish  numerical  values,  from  given  input  information, 

for  V,  (1/V)J,  (1/V)*',  M  and  (T  /T  )  at  that  x*. 

e  we 

3)  Compute  R^  from  Eq.'k)  and  chen  compute  P*  from  (4)  also, 

4)  For  the  given  M  and  (7  /T  ),  ’■'ead  A  from  Figure  1. 

e  •?  e 

5)  Compute  X  from  Ea.(8)  and  hence  compute  X/X 

ra  ^  max 

6)  Enter  Figure  2  or  Table  1  at  (X/X  )  and  read  f*  and  g*. 

7)  Using  A,  f*,  g*  and  R*  from  above,  compute  the  local  rate 
of  change  (dX/dx*)  from  the  basic  relation,  Eq.(7).  You  are  now 
ready  to  "integrate*  forward  to  the  next  position  (x*  +  Ax*)  using 
a  numerical  formula  such  as  X(x*+Ax*)  =  X(x*)  +  Ax*  X’(x*+Ax*). 
Alternately,  the  graphical  "method  of  isoclines"  is  very  convenient, 
since  we  are  dealing  only  with  a  first  order  differential  eauation. 


Numerical  Example: 

At  the  initial  position  of  the  supersonic  relaxing  flow 
experiment  of  F.  V.  Zwarts  (Case  1  in  what  follows),  we  are  given 
the  following  information: 


o 


Cf  =  0.000894 
X  -  0.75  inches 

U  *  U  *  2205  ft/sec 

G  O 

v  =  0.00127  ft2/sec 
e 

V  a  1.0 

(l/V)’  i  0.008  t  .  . 

)  per  inch 

(l/V)"  =  0.003  ~ " 

M  a  4.02 
e 

Adiabatic  Wall 


(L 


1“) 


This  is  sufficient  input  to  proceed  forward  with  Eq.(7),  as  follows: 

Step  1.  Establish  (T  /T^)  from  an  adiabatic  wall  formula  with,  say, 

a  recovery  factor  of  0.89  for  turbulent  flow: 

T  /T  -  1  ♦  0.89(0.2)  M2  a  3.88 

we  e 

Also  estimate  the  viscosity  ratio  from  a  power-law  formula  for  air: 

U/p  i  (T  /T  )0,67  -  2.47 

w  e  we 

Step  2.  Compute  the  initial  value  of  X  from  the  skin  friction: 

X  =  (2/Cf)1/2  =  47.3 

Step  3.  Compute  R^  and  R*  from  Eq.(4):  use  L  =  1  inch  for  convenience. 

R.  a  [2205(l/12)/0.000127]/(3.88)1/2/(2.47)  a  3.56  x  105. 

Li 

R*  a  R^U/V)’  =  1.25  X  108 

Step  4.  For  a  4.02  and  (T  /T  )  =  3.88  (adiabatic),  read  A  =  1.63 
from  Figure  1. 

g 

Step  5.  From  Eq.(8),  compute  X  =  8.7(1.53)  log. ~(1. 25x10  )  a  115. 

TT13X  J.V 

Hence  X/X  a  (47.3)/115.  a  0.412. 
max 

Step  6.  From  Figure  2  or  Table  1,  at  0.412,  read  f*  a  0.294  and 
g*  =  0.011.  For  convenience,  compute  [3g*  =  0.131. 

Stop  7.  Evaluate  the  derivative  from  Eq.(7): 


10 


dx 

dx* 


O^0-13K1-63)21 


0.008 

(1.0) 


{1  ♦  3(0.131)/1.632]  + 

0.16  (0.294)  (1.63) 


3 


0.0091  +  0.121 

0.204 


0.64 


This  is  a  rather  small  rate  of  change ,  corresponding  to  a  modest 
adverse  pressure  gradient.  We  may  integrate  forward  about  three 
inches,  which  gives  about  a  4%  change  in  X.  In  no  case  should  we 
try  to  extrapolate  forward  more  than  about  a  3%-5%  change  in  X  ,  and 
in  this  way  we  will  avoid  serious  numerical  error.  Here  we  estimate: 


X( 3.75")  =  X(0.75")  +  (3.0")  ^*(0.75") 

=  47.3  +  (3.0)(0.64)  =  49.22 

Cf(3.75")  =  2/(49. 22)2  s  0.0008256 

If  Me,  V,  etc.  are  known  at  3.75",  we  are  now  in  a  position  to 
evaluate  (dX/dx*)  at  this  new  station  and  proceed  forward  again. 

If  this  derivative  is  markedly  different  from  the  previous  one  of 
0.64,  we  should  backtrack  and  take  an  "average"  slope  over  the 
previous  three  inches  to  estimate  a  more  accurate  X(3.75").  Note 
that  we  are  using  Ax*  in  inches  here,  since  we  chose  L  *  one  inch. 

If  we  had  chosen  L  =  18"  (the  total  length  of  Zwarts'  model),  then 
both  and  (1/V)'  would  be  eighteen  times  larger  and  hence  R*  would 

g 

still  be  exactly  the  same  (1.25x10  ),  leading  to  the  same  values  of 

2 

f*  and  g*.  However,  (1/V)”  would  be  (18)  larger  and  the  result 
would  be  that  (dX/dx*)  is  eighteen  times  larger.  Meanwhile,  Ax*  is 
eighteen  times  smaller,  so  that  the  value  of  X  at  3.75"would  be 
still  the  same  (49.22).  Thus  the  choice  of  l  does  not  affect  the  result. 
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tv  DIGITAL  COMPUTER  CALCULATION 


The  basic  theory,  Equation  (7),  has  also  been  programmed  for  use 
on  a  digital  computer,  The  suggested  FORTRAN  program  is  on  the  next 
page,  followed  by  a  subroutine  (RUNGE)  which  implements  the  required 
numerical  integration  by  a  Runge-Kutta  procedure. 

The  user  must  provide  a  smooth  curve-fit  to  the  dimensionless 
freestream  velocity  distribution  V(x*).  The  suggestion  here  is  for 
a  fourth  order  polynomial: 

2  3  4 

V  =  a  +  a,x  +  a.  x  +a_x  +  a,  x 
0  12  3  4 


The  first  READ  statement  calls  for  the  user's  coefficients  (a^a^a^, 
a3,a4).  If  some  other  expression  is  used  (exponential,  sine  wave,  etc), 
the  user  must  himself  modify  the  three  statements  following  statement 
ten,  which  compute  V,  VP,  and  VPP.  Immediately  following  these  three 
are  three  more  statements  which  invert  and  compute  (1/V)'',  (l/V)', 
and  (1/V).  These  should  not  be  changed. 

The  second  READ  statement  calls  for  initial  data  appropriate  to 
the  particular  run: 


U  L/v 
o 


e 


AX*,  X* 
o 


X* 

max 


The  program  will  increment  itself  in  steps  of  AX*  (which  is  called  H 

in  the  program)  and  will  stop  when  XMAX  is  reached. 

At  statement  21,  the  user  must  provide  an  analytical  estimate  of 

the  freestream  Mach  number  distribution  M  (x*).  These  need  not  be  a 

© 
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FORTRAN  Program  for  Solving  Ecuation  (7) 


F”  PENSION  Y( 10)  ,  F(  10  ) 

C  READ  IN  POLYNOMIAL  COEFFICIENTS  FOR  V(X*)  AND  ALSO  INITIAL  DATA. 

7  PEAD(5,*)AZ,A1,A2,A3,A4 
READ(5,*)RZEPO,CFZERO,H,X,XXAX 
Y ( 1 )  -  S-:?T(?./CFZERC) 

M  =  0 

8  IF(X-XMAX)  6,6,7 

6  CALL  RfJK~C.,Y,F,X,F,M,K) 

00  TO  (10,20),  K 

C  NOW  COMPUTE  LOCAL  V,  V’,  V”,  PRC ORAM  COMPUTES  1/V,(1/V) ’  ,(1A'J". 

C  BE  CAREFUL,  THOUGH:.  NOTICE  FOR  ZWARTS  DATA  WE  USE  (X-0.75),  NOT  X. 
10  Z  =  X  -  0.75 

V  =  AZ  «•  A1*Z  *  A2*Z**2  +  A3*Z**3  ♦  A4*Z**4 
V?  =  A1  +  2.*A2*Z  ♦  3.*A3*Z**2  •*-  4.*A4*Z**3 
VFP  =  2.*A2  +  6. *A3*Z  +  12.*A4*Z**2 

VFP  =  (2.*VP*VP/V  -  VPP)/V/V 
VP  =  -VP/V/V 

V  =  l./V 

C  NOW  COMPUTE  THE  LOCAL  MACH  NUMBER.  AN  APPROXIMATION  IS  O.K. 
IF(X-6.75)21,22,22 

21  CM  =  4.125  -  X/6 . 

QO  TO  23 

22  CM  =  3.0 

C  NOW  COMPUTE  ADIABATIC  WALL  AND  LOCAL  WALL  TEMPERATURE,  THEN  ,,A‘. 

23  TAW  =  1.  +  0.89*0. 2*CM**2 
7*  =  TAW 

A  =  TW  ♦  TAW  -  2. 

3  =  TAW  -  TV 

C  =  SQRT((TAW  ♦  TW)**2  -  a.*7W) 

A  =  SQF7(?AV-i. )/(ASIN(A/C)  +  ASIN(B/C)) 

C  TEST  TO  SEE  IF  EQ.( 10)SHCULD  BE  USED,  OTHERWISE.  USE  E0-(7). 
IF(VP)41,41,42 

41  F( 1)  =  RL/V/EX?(.48*Y(l)/A)/8.  +  5.5*VP/V 
GO  TO  6 

42  RSTAR  =  RZERC7SQRT(rw)/TW**0. 67/VP 
0  =  Y( l)/8.7/A/AL0G10( RSTAR) 

IF(Q-0.4)41,41,S0 

50  ir(Q-l. )5I,52 ,52 
52  WRITE(6,606)X 

606  FORMATC  SEPARATION  HAS  OCCURRED  AT  X  =  ’,E12.4) 

GO  TO  7 

51  Z  =  1.  -  Q 

GS  =  1.  -  2 . 3*  Z  +  1.76*Z**3 

FS  =  (2.424*2  +  I.443*Z**2)*EXP(-44.*Z**6) 

GR  =  3.*GS*RSTAR**0. 07 

Q  =  V?*(i.  +  3.*SR/A/A)/V  «•  VP?*GR*A*A/V? 

F(l)  =  Q/.16/rS/A**3 
GO  TO  6 

20  CF  =  2./Y(l)**2 

WRITE( 6 ,707)  X.  Y(l),  CF,  V,  CM 
707  F0RMAT(2F10.4,2F14.8,F1D.4) 

GO  TO  8 
END 


1.7 


S U3 ROUTINE  RUNGS  ( N ,  V ,  ?.  X,  H,M.  K) 

PENSION1? (10 )  ^2^0 )  ^ A  CALCULATI0N  BY  GILLB  METRO 

GO  TO  (1X5,3.'’;  ,  M 
1  DO  2  1  =  1.  N 
?  Q  ( I )  =  0 
a  =  c; 

GO  TO  9 
3  A  =  l.^O^lO? 

„  ACCURACY’  USE  A  =  l-'7G'710c'7eii&65-^5?i-ii 

"  EO  6  I  =  1,  n 

-  Y(*)  =  Y(I)  ^  A*(F(I)*K  -  0(1)) 

'  0(1)  =  +  (1.  -  3-*A)*Q(l) 

A  =  .2928932  ' 

“  M0RE  ACCURACY»  SET  A  =  .2928 321 881  ^452^756 

~  ro  8  f=  i.n 

8  «(T0-  V(X)  -  H»F(1)/c.  .  Q(I)/3. 

X  =  2 
GO  TO  10 

9  x  =  1 
1C  RETURN 

ENT 
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smooth  curve-fit.,  since  the  exp^ssion  will  not  be  differentiated  and 
is  used  on^v  to  compute  van  Driest’s  parameter  A  from  Eauation  (5). 

Finally,  iust  below  statement  number  23,  the  local  wall  temperature 
must  he  estimated.  In  the  example  shown,  it  is  simply  set  equal  to 
the  local  adiabatic  wall  temperature  (the  conditions  of  Zwarts*  data). 

Again,  any  reasonable  estimate  will  do,  because  7^  is  used  onlv  to 
compute  van  Driest's  parameter  A,  which  Figure  1  shows  is  onlv  slowlv 
varying  with  both  and  M  . 

All  other  procedures  in  the  theory  are  carried  out  by  the  computer 
program,  including  testing  for  boundary  laver  separation  and  printing 
out  the  separation  point  if  \  exceeds  X 

TTI9X 

Example  1.  THE  SUPERSONIC  RELAXING  FLOW  OF  ZWATTS  (UNPUBLISHED) 

In  unpublished  work.  Dr.  F.  W.  Zwarts  of  McGill  University 

measured  a  two-dimensional  supersonic  flow  which  decelerated  sharply 

from  Hach-4  to  Mach- 3  in  about  eight  inches  and  thereafter  remained 

nearly  constant  at  Mach-3.  This  flow  was  used  as  an  example  to  test 

the  finite  difference  calculations  of  Bradshaw  (B).  It  is  a  difficult 

test  for  any  method.  The  freestream  velocity  V(x)  and  uach  number 

M  (x)  ar'  shown  in  Figure  3.  The  velocity  is  fit  bv  the  solid  line  to 
e 

a  least -scuares  fourth  order  polynonial:  (U  =  2204.  ft/sec) 

o 

V  =  1.0013  -  0.00051  x  -  0.C03885  x2  +  C. 0004466  x3  -  0.0^0013^  x** 

(11) 

The  smoothness  of  the  fit  is  seen  to  be  quite  adeauate,  and  one  can  take 
confidence  in  a  reasonably  smooth  first  and  second  derivative.  The  Mach 
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COMPARISON  OF  THEORY  WITH  THE  SUPERSONIC  RELAXING 
FLOW  EXPERIMENT  OF  ZWARTS  (UNPUBLISHED). 


number  is  fit  adequately  by  two  linear  nieces: 


x  <  6.75":  M  =  4.125  -  x/6 

6  (12) 

x  >  6.75":  M  =  3.0 

e 

Note  that  x  is  in  inches  in  both  Eqs.(ll)  and  (12),  That  is,  L  =  one  inch. 

Also  shown  on  Fieure  3  is  an  exnonential  aoDroximation  to  "  (x): 

e 

V  =  0.9165  +  0.0835  exn[-0.03  x  -  0.05  x2)  (13) 

This  is  plotted  as  a  dotted  line.  The  *it  is  excellent  up  to  x  =  12 
inches,  after  which  the  curvature  in  the  data  -  which  mav  he  an  artifact  - 
is  not  fit  by  the  exponential  curve. 

The  initial  data  to  be  read  into  the  proprar!  are  as  follows: 

R  =  U  L/V  =  ( 2204 )( 1/12 )/(0.C00106)  =  1.73  x  10. 

c  o  e 

CF  =  0.00089“  (*rom  Zwarts'  data). 

‘  o 

x  =  0.75  inches, 

o 

h  -  Ax  =  0.5  inches. 

x  =18  inches, 
m ax 

The  choice  Ax  =  0.5"  is  an  engineering  iudvnent  desimed  to  hreav  the 
interval  0  to  L  down  into  about  forty  steps.  Doubline  this,  to  Ax  =  1", 
will  cause  a  numerical  error  of  about  one  per  cent  in  the  computed  skin 
friction.  The  co."souted  results  fron  the  nropram  are  also  shown  in  Fieure 
3.  Both  the  polynomial  and  exponential  fits  give  theoretical  C^.  In 
reasonable  agreement  with  the  Preston  tube  data.  Both  are  far  more  accurate 
than  the  finite-difference  calculations  of  Bradshaw  and  Femss  (8),  also 
shown,  which  fall  thirty  to  forty  per  cent  low. 
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Example  2.  THE  UPPER  SURFACE  OF  AN  NACA  S4A210  AIRFOIL 


Mann  and  Whitten  (9)  measured  skin  friction  with  a  Stanton  tube  and 
a  Preston  tube  on  the  upper  surface  of  an  J/ACA  64A210  airfoil  in  high 
subsonic  flow,  Ha=  0.7.  The  measured  velocity  distribution  for  zero 
angle  of  attack  is  shown  in  Figure  4.  An  excellent  fit  is  nrovided  by  the 
solid  line,  a  fourth  order  least  squares  polynomial: 

V  =  0.8664  +  2.6494X*  -  5.6028X*2  +  4.1633x*3  -  1.2079x*4  (14) 


Since  the  flow  was  subsonic,  the  static  temperature  was  nearly  constant 
and  thus  the  freestreani  Mach  number  was  nearly  Droportional  to  V: 

M  =  0.70  V  (15) 

e 


For  initial  data,  Mann  and  Whitten  suggest  that  transition  to  turbulence 
occurred  at  about  x*  =  0.1,  with  C^.  somewhere  between  0.003  and  0.004; 


R  =  9.55  x  106 

o 


x*  =  0.1 

o 

h  =  Ax*  =  0.05  (16) 


Cf  =  0.003  or  0.004  ; 

o  x*  =  1.0 

max 

T  i  T 

w  aw 

These  values  were  run  with  the  computer  program  of  rhis  report.  The 
results  are  shown  in  Figure  4.  Both  runs  are  in  reasonable  agreement  with 
the  data.  The  dotted  curve  (Cf0  =  0.003)  shows  separation  at  the  trailing 
edge,  which  was  probably  the  case.  Mann  and  Whitten  themselves  demonstrated 
good  agreement  with  the  Karman  integral  theory  of  Sasman  and  Cresci  (10), 
and  also  with  a  finite-difference  procedure  by  Cebeci  et  al  (11). 
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Example  3.  SUPERSONIC  FLOW  PAST  A  CURVED  COMPRESSION  RAKP 


'll 


Sturek  and  Danberg  (12)  placed  a  faired  curved  raTip  on  the  floor 
of  a  supersonic  wind  tunnel  and  generated  a  smooth  decrease  in  freestreara 
Mach  number  from  3.5  to  2.9  in  a  distance  of  ten  inches.  The  freestream 
velocity  data  are  shown  in  Figure  5.  They  are  fit  quite  well  by  the 
solid  curve,  which  is  a  least  squares  polynomial: 

V  =  1.1132  -  0.34396  x  +  0.005943  x2  -  0.0003263  x3  +  0.00000592  x* 

(17) 

where  x  is  in  inches.  The  derivatives  of  this  fit  may  not  be  reliable, 
since  there  are  only  eight  data  points  to  be  fit  and  none  are  in  the 
important  region  where  the  r»mp  begins.  The  step  size  and  range  are: 

x  *  6  inches,  h  =  Ax  =  0.5  inches,  x  =22  inches.  (18) 

o  ’  max 

Skin  friction  data  are  shown  for  two  different  tunnel  stagnation  pressures, 
which  correspond  to  two  different  Reynolds  numbers  and  initial  values: 

p  -  psia  U  (l")/v  C* 

o  _  o _ e _  *o 

38.7  358,000.  0.001085  ,  .  . 

}  T  *  T 

56.0  538,000.  0.00103  aw 

These  were  run  on  the  computer  program,  and  the  results  are  shown  in 
Figure  5.  It  is  seen  that  the  thecry  is  in  excellent  agreement  with 
the  skin  friction  data.  No  ether  theory  has  beer,  aoplied  to  this  data, 
and  in  fact  most  other  theories,  both  Karmen  integral  and  finite 
difference  type,  are  invalid  for  this  experiment  because  of  the  surface 
curvature,  which  causes  large  pressure  gradients  normal  to  the  wall.  (8.10) 
The  present  thecry  is  insensitive  to  normal  pressure  gradients. 
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Example  FLOW  PAST  A  FLAT  PLATF 


Far  flow  at  zero  incidence  past  a  flat  plate,  a  »  V*  *  V"  *  0, 
and  the  theory  reduces  to  Equation  (20),  which  further  reduces  to: 

~  R.  exp(-0.48  X/A)  ,  with  X  *  X  at  x  *  x  .  (IS) 
dx*  8  L  ^  ’  o  o 

The  variables  are  separable  and  integrable,  and  an  exact  solution  is 
possible: 

1. 

C£  =  0.455  A-2  tn“2[exp(0.48X  /A)  +  %^R  -R  )(|*)2  (^)  1  ,  (20) 

f  -  o  A  x  xq  Tw  uM 

where  R^  =  U^x/u^  -s  the  local  Reynolds  niraber.  This  f rmula  may  be 
used  to  proceed  from  a  position  C^.  (R  )  downstream  to  a  new  value  C^iR^). 

If  the  plate  is  turbulent  from  the  leading  edge,  X  *  R^  *  0,  and  Eq.(20) 
reduces  to  the  flat  date  formula  of  reference  3: 

C,( flat  elate)  *  0.455  A-2  ln~2[^^- R  (T  'T  )1/2(u  /p  )]  (21) 

f  Axewew 

When  compared  in  reference  3  with  almost  all  available  data  on  turbulent 
flat  plate  skin  friction  (657  reported  values),  Eq.(21)  was  shown  to  be 
the  most  accurate  formula  ever  derived  in  terms  of  mean  absolute  error. 

It  is  more  accurate  than  the  tentative  formula  we  first  suggested  (in  ref.  2) 
and  is  a  consequence  cf  a  better  correlation  of  the  integral  functions  G  and 
H  of  that  reference. 

Some  numerical  values  from  Eq.(21)  for  adiabatic  walls  are  given  in 
Table  2.  and  a  plot  showing  the  effect  of  wall  temperature  is  given  in 
Figure  6.  Note  that  there  is  a  substantial  effect  of  Reynolds  number  in 
this  plot,  so  that  the  common  practice  in  the  literature  of  plotting  all 
data  on  a  single  figure  of  this  type  can  be  very  misleading. 


Table  2 


TURBULENT  SKIM  FRICTION  OK  ADIABATIC  FLAT  PLATES,  EQ.  (21} 


LOCAL  REYNOLDS  HUMBER.  F 


MACH 

HUMBER 

xJ 

1C6 

9 

0.006012 

0.003759 

1 

0.005703 

0.003522 

2 

0.005049 

0.003025 

4 

0.003892 

0.002146 

6 

0.003225 

0.001620 

8 

0.C02868 

0.001308 

1C 

0.002696 

0.001110 

x 


ID7 

1C8 

• 

10  ^ 

0.002570 

0.001866 

0.001419 

0.002390 

0.001727 

0.0013HP 

0.002012 

0.001435 

0.001074 

0.001357 

0.000935 

0.000683 

0.0C0972 

0.000647 

0.000462 

0.000745 

0.0C1491 

0.000335 

C. 000601 

0.000376 

0.000259 

Three  other  examples  of  the  use  of  this  new  tneory  to  predict 
turbulent  skin  friction,  particularly  in  supersoni flows, are  given  in 
references  3  and  E . 


V  CONCLUSION'S. 

It  has  been  shown  that  one  can  compute  turbulent  skin  friction  in 
two-dimensional  flows  with  arbitrary  conditions,  using  Eo.(7)  of  this 
report,  which  has  been  programmed  for  ready  use  with  a  digital  computer. 

The  user  must  only  supply;  1)  an  initial  value  of  skin  friction;  2)  a 
smooth  curve  fit  to  the  freestream  velocity  distribution;  and  3}  reasonable 
estimates  of  the  freestream  .'4ach  number  and  wall  temoerature  distributions. 
The  authors  believe  this  theory  is  the  rnost  accurate  in  the  literature. 
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